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CN , We examine the constraints imposed by the requirement of successful nucle- 

I osynthesis on models with one large extra hidden space dimension and a single 

^ I bulk neutrino. We first use naive out of equilibrium conditions to constrain 

the size of the extra dimensions and the mixing between the active and the 
bulk neutrino. We then use the solution of the Boltzman kinetic equation 
. for the thermal distribution of the Kaluza-Klein modes and evaluate their 

contribution to the energy density at the BBN epoch to constrain the same 
parameters. 



I. INTRODUCTION 



There has been a great deal of interest and activity in the last two years on the possibility 
that there may be one or more extra space dimensions in nature which have sizes of order 
^ ■ of a millimeter This has been driven by the realization that string theories provide a 
' completely new way to view a multidimensional space time in terms of a brane-bulk picture, 
where a brane is lower dimensional space-time manifold that contains known matter and 
forces and the bulk consists of the brane plus the rest of space dimensions where only 
gravity is present. The resulting picture replaces Planck scale by the string scale as the new 
fundamental scale beyond the standard model. The relation between the familiar Planck 
scale and the string scale is given by the formula 

Mj., = Ml+^R.R^R^ ■■■Rr, (1) 

For i?i ~ i?2 = -R and i?3 ~ i?4 ~ ■■ = M~^, this relation leads to i? ~ millimeter for 
M^, ~TeV. The fact that the familiar inverse square law of gravity allows for the existence of 



*e-mail: hsgoh@glue.umd.edu 

^ e-mail : rmohap at ©physics . umd . edu 



such sub-millimeter size extra dimensions made these models interesting for phenomenology 
0. An added attraction was the fact that a whole new set of particles are present at the 
TeV scale making such theories accessible to collider tests. Furthermore since there are no 
high scales in the theory, there is no hierarchy problem between the weak and the Planck 
scale; this provided an alternative approach (resolution ?) to the familiar gauge hierarchy 
problem. Obviously, the picture would become much more interesting if collider experiments 
such as those planned at LHC or Tevatron fail to reveal any evidence for supersymmetry. 

Even though these models present an attractive alternative to the standard grand uni- 
fication scenarios, there are two arenas where the simplest TeV string scale, large extra 
dimension models lead to problems: (i) one has to do with understanding neutrino masses 
and (ii) second is in the domain of cosmology and astrophysics. 

The reason for the first is that the smallness of neutrino masses is generally thought to 
be understood via the seesaw mechanism the fundamental requirement of which is the 
existence of a scale > 10^^ or 10^^ GeV, if the neutrino masses are in the sub eV range. 
Clearly this is a much higher scale than of the TeV scale models. A second problem 
is that if one considers only the standard model group on the brane, operators such as 
LHLH/M^, could be induced by string theory in the low energy effective Lagrangian. For 
TeV scale strings this would obviously lead to unacceptable neutrino masses. 

In the domain of cosmology, the problems are related to the existence of the Kaluza- 
Klein (KK) tower of gravitons generally which lead to overclosure of the Universe unless the 
highest temperature of the universe is about an MeV 0]. This can not only cause potential 
problems with big bang nucleosynthesis but also with understanding of the origin of matter, 
inflation etc. There are also arguments based on SN1987A observations that require that 
M, > 50- 100 TeV §. 

The neutrino mass problem was realized early on and a simple solution was proposed 
in ref. 0. The suggestion was to postulate the existence of one or more gauge singlet 
neutrinos, z/^, in the bulk which couple to the lepton doublets in the brane. We will call this 
the bulk neutrino. After electroweak symmetry breaking, this coupling can lead to neutrino 
Dirac masses, which are of order hv^k^*/ ^pt- This leads to niy ~ 10~^ eV. The dominant 
nonrenormalizable terms have to be forbidden in this model. The simple way to accomplish 
this is to assume the existence of a global B-L symmetry in the theory. The only difficulty 
with this assumption is that string theories are not supposed to have any global symmetries 
and one has to find a way to generate an effective B-L symmetry at low energies without 
putting it in at the beginning. 

There is an alternative scenario for neutrino masses where one abandons the TeV 
string scale but maintains one large extra dimension and avoids the problem associated with 
nonrenormalizable operators. The relation in Eq. (1) then gets modified to the form 

Ml, = MlR (2) 

Since the string scale in these models is in the intermediate range i.e. 10^ GeV or so, the 
cosmological overclosure problem is avoided. The active neutrino masses in such models 
could arise from seesaw mechanism or from the presence of bulk neutrinos. The inclusion of 
the bulk neutrino however brings in new neutrinos into the theory which can be ultralight 
(i.e. 1 where R is the size of the large extra dimension) and can play the role of the 
sterile neutrino, which may be required e.g. if the LSND results are confirmed. 
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Both the above approaches have the common feature that they introduce a bulk neutrino 
into the theory, which is equivalent in the brane to an infinite tower of sterile neutrinos. All 
of these neutrino modes mix with the active neutrinos in the process of mass generation 
PJ^. For extra dimension size of order ~ mm, the KK modes have masses typically of order 
nR~^ ~ n X 10~^ eV or so, where n = 0, 1, 2, 3, ■ ■ ■. The presence of this dense tower of extra 
sterile states coupled to the known neutrinos leads to a variety of new effects in the domain of 
particle physics and cosmology |]TU|-[T^, which in turn impose constraints on the allowed 
size of the extra dimensions. In this paper, we focus on the cosmological constraints that 
may arise from the contribution of the neutrino states to big bang nucleosynthesis. 

The constraints from big bang nucleosynthesis on bulk neutrinos were considered in 



Ref. where the cases of > 2 extra dimensions were discussed. As it was noted 

there, light KK modes of the bulk neutrinos could easily be produced in the early universe, 
when temperature is of order an MeV or more. Thus, there is the danger that they could 
make large contribution to the energy density of the universe at the epoch of big bang 
nucleosynthesis (BBN) and completely destroy our current, successful understanding of the 
primordial abundance of He^, D and Li^ |T^. In particular, present abundance data from 
metal poor stars is well understood provided we do not have more than one extra active 
neutrino in the theory in addition to the three known ones i.e. {i'e,i'^,i'r)- By the same 
token, if there are extra species of neutrinos that do not have conventional weak interactions, 
their masses and mixings to known neutrinos must obey severe constraints [|14| . 

Our goal in this paper is to revisit this issue in the context of models with only one large 
extra dimension. The first reason for undertaking this analysis is that the class of models 
with intermediate string scale ~ 10^ Gev and one large extra dimension |]^,|^ have certain 
theoretical advantages and they are also free of the cosmological and astrophysical problems 
that seem to plague the TeV scale models. Secondly, the number of KK modes in this case 
are much fewer than models with larger number of large extra dimensions and therefore, 
one would expect the constraints to be somewhat less restrictive. 

We also wish to emphasize that one large extra dimension could also occur in models with 
string scale in the 100 TeV range, where one can satisfy the Planck-scale-string-scale relation 
in Eq. (1), if the compactification is not isotropic, e.g. for a string scale of 100 TeV, if two 
extra dimensions have sizes r ~ GeV~^ and one has R ~ millimeter, i.e. Mj,^ = M^Rr^. 

We organize this paper as follows: in sec II, we discuss the constraints of big bang 
nucleosynthesis on the size of the extra dimension in the presence of the bulk neutrino and 
the mixing of the bulk neutrino to the active one using simple out of equilibrium condition. 
In section III, we use Boltzman equation to study the generation of the bulk neutrinos from 
active neutrino interactions in the early universe and find the constraint of BBN on the same 
parameters as in sec. II. The numerical calculations leading to our final results are given in 
sec. IV. In appendix A, we explain the details of the out of equilibrium condition for the KK 
modes of the neutrino. 



II. CONSTRAINTS FROM OUT OF EQUILIBRIUM CONDITION 

The class of models, we will be interested in, are assumed to have one large extra di- 
mension with a single bulk neutrino ub, which means that masses of the KK modes of ub 
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are integer multiples of the basic scale /i = ~ 10~^ eV. This is one of the parameters 
that we expect the BBN discussion to constrain. The second parameter is the mixing of the 
KK modes with the active brane neutrinos, e.g. z/g. It is true in both classes of models i.e. 
both TeV and intermediate scale type, that the typical mixing parameter of the active 
neutrino to the nth KK mode scales like 9en — the range of interest for phenomenol- 
ogy. The parameter 6 depends on the size of the extra dimension and other parameters 
of the theory such as the weak scale etc. For instance, in TeV scale models one has 
6 ~ V2hv kM^R ^i^greas in the local B-L models, the relation is 6^ ~ hih^mR where we have 
chosen ~ vr. In general, therefore, BBN discussion will give a correlated constraint 
between /i and 9. Obviously, for ^ = 0, there is no BBN constraint on or the size of the 
extra dimension. 

We will also assume that the universe starts its "big bang journey" somewhere around a 
GeV or so and when it starts, the universe is essentially swept clean of the sterile neutrino 
modes. This can happen in inflation models with a low reheat temperature. We choose such 
a low reheat temperature essentially for reasons that in models with large extra dimensions 
higher temperatures would lead to closure due to production of graviton KK modes. 

To see the origin of constraints, let us note that at high temperatures (i.e. T ^ MeV's), 
there are two ways the KK modes of the sterile bulk neutrino can be created: (i) first, 
neutrino scattering and annihilations and (ii) the oscillation of the active neutrinos into the 
sterile KK modes. It is important to stress that in building up the oscillation, the scattering 
process is important, since otherwise there will be back-and-forth oscillation and no build-up 
of the sterile modes. 

Since there is an infinite KK tower of these neutrino modes, the higher the temperature 
the larger the number of modes that can get created. Once these modes are created, they 
may decay or annihilate to produce the lighter particles (lighter neutrino modes or KK 
modes of the graviton etc). In general, it is reasonable to expect that this process of decay 
or annihilation will not be efficient enough 0] to eliminate all the KK modes. As a result, 
many of them will stay around at the BBN temperature and contribute to the energy density. 
The present understanding of the big bang nucleosynthesis [|l^] relies on the assumption that 
the total energy density at the BBN era is pbbn = fofl'*^'^ with g* = 10.75 coming from 
the contribution of photon, e+e", and the three species of neutrinos. The uncertainties in 
our knowledge of the He^, D2 and Li'^ content of the universe allow that one could have 
g* ~ 12.5 (or one extra species of neutrino). We will require that any additional contribution 
to Pbbn coming from the bulk neutrinos generated at higher temperature be less than the 
contribution to Pbbn equivalent to one extra species of neutrino. 

The first step to ensure this would be to enforce the condition that the production rate 
for any singer bulk neutrino mode is less than the Hubble expansion rate of the universe 
between 1 GeV to MeV. This will prevent any bulk neutrino mode from equilibrium. The 
equilibration of any neutrino mode is unacceptable because it can contribute more than 
allowed energy density to the universe. Clearly this condition cannot be a completely reliable 
method of constraining the parameter space when there are such a large number of final 
states present since a small oscillation into each mode can ruin BBN results. However, we 
use this as a "warm up" to our final (hopefully more refined treatment) and as a basic 
guideline for more precise and stronger condition. We present the details of this discussion 
and its results in this section fully realizing that this is not going to be the "final story" . Our 
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next step to obtain constraints is to use the Boltzman equations for the time evolution of the 
density of the KK modes and obtain their contributions to the energy density at the BBN 
era and demand that this energy density is less than that of one extra species of neutrino. 
We will discuss this in the next section. We find it intriguing that both these methods lead 
to bounds which are very close to each other. 

We begin with some pedagogical comments by considering one sterile mode which mixes 
with one of the active neutrinos and by ignoring the matter effects. We present this discussion 
as a way to appreciate the importance of the matter effects. In this simple case, it turns out 
that it is the oscillations rather than direct production in scattering that give the strongest 



bound |T5| 



To see that direct production is not so important, note that the the production rate 
for the nth KK mode is roughly given by G^pT^ ^^9"^ . Note that dependence on the mode 
number cancels out. Setting this < the Hubble expansion rate , we get, for the decoupling 
temperature T* ~ g*^/'^ /{Gpfi'^O'^Mpi) ^ 10^^^^. Thus below this temperature all produc- 
tion processes for arbitrary modes are far out of equilibrium. Thus we only have to consider 
ub production through oscillations. 

The transition rate to one sterile (or KK) mode from oscillation in the absence of matter 
effect can be estimated as follows: the oscillation is a quantum mechanical phenomenon that 
gets interrupted as soon as a collision takes place. The amount of time which allows a build 
up of the oscillation to the sterile state is therefore given by the time between the collisions 
T i.e. inverse of the collision rate, W^k — G\T^ . The probability for transition to a sterile 
mode in time r ~ W~l is given by 

P ~ sin^2een sin^ {—\ (3) 



osc 



where oscillation time Tosc — ir-^ — —3 — • If time between collisions of the z/e, denoted 
by Tyjk, is much less than the oscillation time, this expression is simplified and we get the 
production rate for the sterile KK modes 

Wkk ^ PW^k ^ stn'2e,n (^) (4) 

V '^osc / 

This condition i.e. r^k = Gp^T~^ < E/Am'^ is satisfied around 2-3 MeV for m < 10~^ eV. 
Below this temperature we must approximate the oscillation factor by 1/2. Staying above 
3 MeV and ignoring the matter effects, we can implement the BBN constraint by requiring 
that the production rate for a KK mode be less than the Hubble 
epoch: 

Wkk < g*'"^ (5) 



If a given KK mode is much lighter than an MeV, we get the bound in literature |T5[ that 

<KOln < 10"^' eV' (6) 

Taking matter effects into account essentially amounts to replacing the vacuum mixing 
angle by the matter mixing angle given by 
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(7) 



where z = -6.3T ^^/"f ^ _M55xioli(T^ ^ f^_^^ ^^^^ ^ 



have assumed that there is no lepton asymmetry. Clearly a{T) is a function of temperature. 
Note that z is always negative. This means that there is no resonance type behavior for the 
mixing angleQ, which makes the calculation more reliable. 

The out of equilibrium condition in the presence of matter effect is given by: 

10^'"^ ^ 2.5V5^r-3io- (8) 

where all mass parameters are in units of GeV; /i = is the KK scale, discussed before. 
This condition must be satisfied for all k and within our temperature window of one MeV 
to one GeV. We include the detail analysis of this part in Appendix A. For 9, /i, T fixed, the 
equation has to be satisfied when the expression in the left hand side of the above expression 



is a maximum. This occurs when k = ya{T) for a{T) > 1. For a{T) < 1, the maximum 
occurs when k = 1. In the first case, we get the limit, 

/i2W2^ < 2.15(^^)^Vrl02e\/2 (9) 

The best bound arises when T is minimum consistent with the condition a{T) > 1 or 
equivalently 

(^)a > AttXL. (10) 

^GeV^ - 2.15 X IQii ^ ^ 

In the second case, we get the limit, 

sin^29 <2.5{-^)-'^^/^{l + a)h0-'^ (11) 

The best bound arises when T is maximum consistent with the condition a{T) < 1 or 
equivalently 

. r ^« < (12) 



^GeV 2.15x1011' 

The temperature should also be within our chosen window of one MeV to GeV. This gives rise 

to three cases. To see the various cases, let us define the temperature Ti = (^ -jiixio'^^ ) GeV 
. The various cases then correspond to (i) Ti < MeV; (ii) GeV < Ti; (iii) MeV < Ti < GeV. 
The limits for various cases are found to be as follows: 



Note that if there was a significant initial lepton asymmetry, the situation would have been 
very different and whether there is a resonance would depend on whether the initial particle is a 
neutrino or an antineutrino. 
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(i) In this case fi"^ lies in the range //^ < 2.15 x 10 ^eV^ 
and the constraint on 3171^29 is given by 

/i^ sin^2e < 7.07 x lO'^el/l (13) 

As sin'^26 is less than 1, this imposes no further constraint on our parameter space. 

(ii) For this > 2.15 X 10^^ eV^ Which is far beyond our range of interest. 

(iii) /i^ is in the range 2.15 x lO^^eV^ < ^2 ^ 2.15 x 10" eV^ 
and the bound is 

fx sin^2e < 5.804 x 10~^e1/ (14) 

The bound of the allowed parameter space due to the out of equilibrium condition is 
adequately described by eq. (p!4D within the range of interest. We plot the result in Fig. 1. 
It corresponds to the small dashed line, the right side of which is forbidden. 



III. BOLTZMAN EQUATION AND BBN CONSTRAINTS 

In this section, we employ the Boltzman equations to get the constraints on /i and 9. Our 
procedure is to calculate the distribution function for the sterile KK modes produced in the 
matter oscillation of i^g, including any possible depletion of their density due to decays all 
the way down to the BBN temperature. We then calculate their cumulative contribution to 
the energy density p at the BBN epoch and demand that this be less than the corresponding 
contribution of one extra species of neutrino. This is the procedure followed in |T7| , p!T|| . 

To estimate this new contribution to energy density pbbn, we have to calculate the 
number of KK states produced at a given temperature. Let us denote by fk = fkiPyt), the 
distribution of the kth mode of the neutrino at the epoch t. The KK modes are not in 
equilibrium at any epoch. The time evolution of the nonequilibrium density is governed by 
the Boltzmann equation given below [|Tl|] [|1^ . 



Ot op J^k Tk 

where we have neglected the contribution of the pair annihilation of the KK modes in the 



right hand side, since it is a very small effect. In Eq. (|T^), t (time), p (momentum) are 
the two independent variables with k, p and sin^26 fixed. Ek{p) = y/p'^ + m^^. H is the 
instantaneous Hubble expansion rate and is clearly a function of time H{t). r(z/Q — > 
^k) = Tfc, the production rate of bulk neutrino, is given by (r/2) (P(z/q — > u^)) where 
(r/2) = 2Gp^T^ is half the interaction rate of the active neutrino in the thermal bath. 
Taking matter effects into account, the probability P is given by 

^^^'^^ ^ ""'^^ " 2 1-2.cos2g. + .2 - (16) 

Note that we have used the averaged probability P to get rid of the momentum dependence 
of the production rate. We will only use this explicit form in our numerical calculation. 
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The averaging is not necessary for the general solution we are going to find later in this 
section. In the above equation, f^^ is thermal distribution of active neutrino i>a- is the 
lifetime of the k-th mode in its rest frame, and = total decay width. For small k, the 
dominant contributions come from the partial width of uj. — >■ decay which is given by 
sin'^ekGjml/ 19271^ = sin^2ekG)^r'k^ / {A x 1927r3) |ni. For big k, it is from YX~=i^k 



Vk'hk-k' ~ — \) /\2tx ~ fc^/i^/127rMp;. We included both contributions in our total 
decay rate. To make the solution more general, we now let all of the functions in Eq. (|T3p 
except H{t) depend on both p and t. 

The general solution of eq. (^) without knowing the form of functions introduced in the 
equation and Ck^i = is found to be 

fkM = fVk{v\x,t),x)USv\x,t),x)e-^'^<^'^dx (17) 
Jti 

provided 

p'(a;,t) =pe/>("')''"' (18) 
Note from Eq. (18) that p'{t,t) = p. 

f* 1 
a{x,t) = / =dx' (19) 

Tk{p{x',t),x')^p^{x',t)+mk^ 
The above solution can be checked easily by the observation that 

(|-^fp|)p'(x,t) = (20) 

and 

(|-^p|W(^,^))=0 (21) 

for any function F{p'p{x,t)) with no explicit dependence on t and p. The time derivative 
on the integration upper limit gives the first term on the right hand side of eq. (|T5|). While 
acting on the upper limit of a{x,t), it gives the second term on the right hand side of eq. 
(|r^). The total energy density of bulk neutrino is then given by 

poo roo fp, r)2 

p{t) = I I -^fk{p,t)dkdp (22) 

in the continuous k approximation. In the next subsection, we discuss the numerical results 
that follow from the above discussion. 



IV. NUMERICAL RESULTS 

In this section, we calculate the total energy density of the universe at T = IMeV from 
(p^). We use temperature (T) instead of time in the integration and proceed in two different 
ways. 
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In the first method, we make some approximations to simphfy the integral in Eq. (17) 
and calculate analytically as far as possible and then numerically estimate the final integrals 
that follow from it. 

In the second method, we evaluate the entire integral in Eq. (17) numerically . The two 
results agree approximately with each other. 

In the calculation, we use H{T) = 1.66^/^^- and t = 0.301(7* ^ as given by Standard 



cosmology. After introducing dimensionless variables p' = f^, x = -^j and a = A;^,and using 
the notation p instead of p' in the final form, Eq. (p2D reduces to 



p(Tf) =A ' / — ^-^^^=^^e''J^dpdadx, (23) 



where 



/^(l + z^-2z./l- 



A = ^:^i^stn^2eGlMpiT^; (24) 



a 



2 2^2 



D = ^ ^a^ + ^^a^]K(x,p) 

y/Wfb^ 4xl927r2 127rM„V ^ ^'^^ 



(25) 



where i^(x,p) = {VTTp^ - p^ + ^ - pHn ^'+v^l^ =). 



"3 , Wp2 + i^^ 



^ = 1 ■ (26) 



and c = j^^yj, ^ = /^/^/- order to simplify the calculation, we also make the following 
approximations. 

• We assume that the effective degree of freedom g^, is not affected by the density of the 
bulk neutrino during the production of bulk neutrino. is given by the standard model 
of cosmology and approximated by a step function as follows: (7* = 61.75 from T = IGeV 
to 200 MeV, = 17.5 from T = 200 MeV to 100 MeV, and g^ = 10.75 from T = 100 
MeV to 1 MeV. Here we actually approximate the degree of freedom as a step function of 
temperature. 

• We use 5m'^ = ml — ml ~ /c^/i^ 

Our Sm'^ are always positive and so lepton number generated from neutrino oscillation 
with nonzero initial lepton number is small [|1^]. We can ignore the contribution of lepton 



number to the Ue potential in matter. This gives 

^6 
J2 



-0.215589 X 10"^^^ (27) 
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c is a constant of 0(1) in the step function approximation of g^,. 

In order to use the analytical method we make some further approximations. Note that 
the integration in Eq. (23) is suppressed by the two exponential functions: one from the 
decay term {D term) and a second from the term. We can therefore cut off the 
integration when either one of them gets smaller than e~^^^. We want to extract some of 
the regions of a and p space in which the J-' term decays much faster than the D term 
within the parameter space of interest and so the latter can be ignored. However, it is 
easier to do it another way. As the contribution from the region where < e~^°° can be 
neglected, we have to examine only region where pa < 100 to make sure D term can be 
approximated to 1. This is same as finding the region where pa < 100 and both ^ x 10^^^ 
and 64fc^/i^ sin^ 29 x 10^^^ less than 0.1. It is easy to see that the only region that may not 
satisfy the above condition is 

. ..AilO^^i , 10^8 ,1, , 
a > rmnli^)-., ( g,^, ,g )^] - (28) 

and 

p < < 1 (29) 

Now we can calculate the total energy density by splitting the integration into three 
parts : 

(i) p > 1 

(ii) p < 1 and a < kmin 

(iii) p < 1 and a > kmin 

We will include the decay term only for (iii). We also simplify our analysis by substituting 
JF = |e~"*' with a possible error of factor 2. For simplicity, we will now treat (7* as a constant 
and so c = 1. Setting cos26 = 1 affects only the term due to the matter effect. This is a very 
good approximation as long as z is negative. For cases (i) and (ii), Eq. ( [23| ) now becomes 



P=^i I i e--'dpdadx (30) 



The X integration can now be done analytically, we approximate \/l + p'^ with p for 
case (i) where p > 1 and with 1 for case (ii) where p < 1. With this simplification, the p 
integration can also be done analytically. For a, we can always divide it into three parts: 
a > 10 , 0.01 < a < 10 and a < 0.01 since kmm > 10^. we will leave the factor to 
the end of the discussion. Now we only look at the integration. The results was summarized 
below 

(i) p > 1 . 

For a > 10, we have 

10^ /-oo 

/ e-'^kdk = 83233 (31) 

6 Jio 

For 0.01 < a < 10, we have 
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flO lilt 1 tan~^ 0-465 



For a < 0.01 



i/ 



0-01 109 TT 1 tan-if 

e 0:216 + 2fc^ ~ ~ ^ ' 

+ — In ^ ^ - (34) 



0.432 4a e 2a 

1 /-o.oi tan-i ^ 



2a Je k 



dk (35) 



(ii) p < 1 and a < k^, 
For a > 10, we have 



For 0.01 < q; < 10, we have 



/ -dk = In 36) 

3 Jio A; 3 10 ^ ^ 



For a < 0.01 we have 

— / + 7^^ o -dk (38 

36 0.216 2a A;2 + a2 a ^ ^ 

This is a very small number which can be neglected. 

(iii) In this region, the decay term have to be included. However, we can remove the p 
dependent by setting p = 1 in the function K{x,p), defined above. The matter effect can 
also be neglected as kmin > 500. After the p integration, with the fact that e"'^™™ < 1 we 
have the upper bound 

/lO"^ POO 

In the above equation, a and f3 are coefficients dependent on kmin- One of them is 0.1 
and the other < 0.1 by the definition of kmin- We also use k = kmin in the result of the p 
integration. For the upper bound, we can just take either the a or the (3 term whichever 
have the value of 0.1. this gives 

^m^^Ki^^ (40) 

4(0.1i^(a;,l))3 



or 



,.o-^,r(i.0.1A-(..l))^ (41) 
Ji 3{0.1K{x,l))-3 
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both give numerically order of 10^. Our result from other parts give about 10^°. This part 
contributes only a few percent of the total. We can use the upper bound 4.6 x 10*. 

The total energy density of bulk neutrinos is the sum of all above , which is found to be 

A , g vr 0.01 f-tan-i^ 1 /-O-Oi tan"^ f „ 10^ kmin, , s 

P = ^ 8.785 X 10^ + — In ^ ^ - — / -^dk + — In ^ 42 

^ 4a e 2a 2a Je k 3 10 ^ ^ ^ 

Compared to the z/g equilibrium energy density at T = IMeV which is 2.824 x lO^^^(G'eV^)'^, 
the constraint that effective neutrino degree of freedom should be less than one, gives the 
constraint on the parameter space 

u \ 

-V- sin^ 29 < 2.5 x 10"^ (43) 

eV J 



We also numerically integrate e.q ([2^ ) without making any simplification except the step 
function g^, and the positive 5m^. Both results from the Boltzman equation and the out of 
equilibrium condition are shown in Fig.Q.The extra effective degree of freedom equal to 1 is 
shown in Fig. 1 as a solid line. The numerical fit is obtained to be 



^.0.92 

eV 



sin^2e < 7.06 x 10"^ (44) 



for /X < leV . 



Log [Sin^2e] 




FIG. 1. The solid and the long dash line are the numerical results which represent the effective 
degree of freedom equal to 1 and 0.1 respectively. The small dash line separate the equilibrium 
and out of equilibrium region. The parameter space below the lines are the allowed region. 



We note that the bounds derived in both ways i.e. out of equilibrium condition and 
Boltzman equation are very similar. 

In Fig. 2, we give the contributions of individual modes to the total energy density as a 
function of the mode number for sin^29 = 10~^ and for a range of values for = fi from 
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GeV^ 




Log [k ] 



FIG. 2. The various lines denote the contribution of the individual modes to the total energy 
density as a function of the KK mode number for various values of inverse size (in eV) of the extra 
dimension. The right-most curve corresponds to fi = 10~® eV and for each successive curve, ^ goes 
up by a factor of 100. 



10 ^ eV to 10^ eV in increasing steps of 10^ eV. The total contribution is the integral over 



each line for a given /i. 



V. CONCLUSION 



In this paper, we have studied the constraints of big bang nucleosynthesis on the models 
with one large extra dimensions and a single bulk neutrino. We find that these bounds 
allow a range of radius of extra dimension and mixing of bulk to active neutrinos that is 
of interest for studying solar neutrino oscillations f^. Our results complement the work of 



Abazajian, Fuller and Patel |T1[] who have derived bounds for the case four, five and six 
extra dimensions. We find it intriguing that the bounds obtained from the naive and crude 
"out of equilibrium" conditions (given in Eq. (14)) are so similar to the ones obtained from 
more detailed considerations based on the density matrix equations. 
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APPENDIX A 



In this appendix, we give details of the out of equihbrium condition used in section II to 
get the bounds on /x^ and mixing 6. Let us first recall that the out of equilibrium condition 
is given by: 

10^^^^ ^ 2.5V5^T-10- (45) 

where all masses are in GeV units and 

2.15 X lO-^T^ a 



Putting in the numerical value of Mpi and we get 

sin'^2e . ,,66.07 



(46) 



sin\^—k') < 2.5v^r-'^10-'" (47) 



(A; + f)2 ^ a 

This condition must be satisfied for all k and within our temperature window of an MeV to 
one GeV. For 9^ /j, T fixed, if the inequality is satisfied when the expression in the left hand 
side is a maximum, clearly it is then always satisfied. The k dependent part is given by 

szn'C-^k^) ^ .m^(66.07y) 

y ^ - (49) 

a 

Here a is independent on k. In this expression, S2n^(66.07y) oscillates so fast that in cal- 
culating the maximum, we can set sin^(66.07|/) = 1 and finding the maximum of ^^j^rjg- 



The maximum is found at y = 1. y = 1 implies k = ^Ja. Because k must be positive 
integer, wc should choose k equal to the closest integer of ^/a. if a > 1. For a < 1, since 
y > Ic and it is easy to see that the maximum occurs at the minimum of y, /c = 1 is the 
maximum point. It is worth noting that, the condition have to be satisfied within a range 
of temperature (and so a continuous range of a). Although at some temperature with a well 
chosen /i, the constraint can be weaker than those we use, small range of temperature of 
MeV will cover the period of sm^(66.07y). So the constraint will be that obtained just by 
setting sin^(66.072/) = 1. 

Now we discuss the two different cases separately, 
(i) For q; > 1 ^ T > ( 2 15x10-7 )^ = 7i-we have constraint 

ii^ sinH9 < 2.15{-^f^/^10'^eV^ (50) 
GeV 

2 1 

(n) For a < 1 ^ T < ( 2 15x10-7 )^- We have constraint 

sin^29 < 2.5(1 + af{-^)-''^^/^10-^^ ~ 2.5{-^)-^^/^lQ-^^ (51) 

GeV GeV 
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In case 2, we have approximate a ^ 1. Without making this approximation, the result 
will give approximately a factor of 2 to sin'^26. Under this approximation, three different 
situation appear, they are (i) Ti < MeV /x^ < 2.15 x IQ-^eV^ ,(ii) MeV < Ti < GeV =^ 
2.15 X 10-^e1/2 < ^2 ^ 2.15 x 10^^ eV^(iii) Ti > GeV ^ /x^ > 2.15 x eV\ 

(i) Ti < MeV 

We have to use eq. (|50|) with T = IMeV and g = 10.75. This will give 

^?sin^2e < 7.07 x lO^^el/^ (52) 

(ii) MeF < Ti < GeV 

We have to use eq.® for T e (Ti, Gel/) and use eq.® for T G (MeV, Ti) with T = Ti 
for both cases. This gives respectively 

/i sin^2e < 2.322 x IQ-^eV (53) 

and 

/i sin'^2e < 5.804 x lO'^eV (54) 

It is easy to see that even without the approximation that a < 1, the second equation can go 
up as high as 1.78 x lO^^el^ when Ti > 1.2MeV. It is still more stringent than the first. We 
use the second equation with ^ = 5. This will at most give another factor (1.5)^ which 
depend on /i. 

(iii) Ti > GeV 

Use eq.(|5TD with T = IGeV and gives 

sin^29 < 1.96453 x 10"^ (55) 

Note that we have omitted the change of the Sm^ due to the matter effect which will only 
increase the as our z is always negative. Increasing will increase the frequency the 
the oscillation term sin^ (66.07?/) and make it easier to be set to 1. 
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